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CRT was a widely used algorithm in the development of watermarking 
methods. The algorithm produced good image quality but it had low 
robustness against compression and filtering. This paper proposed a new 
watermarking scheme through dual band IWT to improve the robustness and 
preserving the image quality. The high frequency sub band was used to index 
the embedding location on the low frequency sub band. In robustness test, 
the CRT method resulted average NC value of 0.7129, 0.4846, and 0.6768 
while the proposed method had higher NC value of 0.7902, 0.7473, and 
0.8163 in corresponding Gaussian filter, JPEG, and JPEG2000 compression 
test. Meanwhile the both CRT and proposed method had similar average 
SSIM value of 0.9979 and 0.9960 respectively in term of image quality. The 
result showed that the proposed method was able to improve the robustness 
and maintaining the image quality. 
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1. INTRODUCTION 

The increasing number of digital media violation has increased the concern about the digital media 
protection. One of the ways to resolve the issue is the watermarking method which becomes rapidly 
developed in the last decade. Watermarking is a way of protecting the digital media by inserting certain 
information called watermark into the media [1]. The watermark contains copyright information that can be 
used as proof of legal ownership. Watermarking methods can be implemented on most of digital media such 
documents, audio, image, and video. In current development, image becomes the most used popular medium 
in watermarking due to its capability of keeping the watermark invisible to the human eye especially in 
heterogeneous area. 

Watermarking methods can be classified into spatial and transform domain methods. Spatial domain 
method is mean of embedding the watermark on the image pixels directly and it has lower complexity but 
less in robustness. Alternatively, the transform domain methods are embeds the watermark on the 
transformed pixels called coefficient in various frequency levels. They have better robustness for retaining 
the watermark [2]. Several development of watermarking in transform domain methods have been conducted 
such Discrete Cosine Transform (DCT) [3]-[5]. Discrete Wavelet Transform (DWT) [6]-[9], Integer Wavelet 
Transform (IWT) [10]. 

DCT is a basic transform method that widely used but it can degrade the image significantly when 
embedding large amount of watermark. The more advanced method DWT is able to represent an image into 
visual feature in multiple level of resolution which can produce watermarked image with better 
imperceptibility and robustness [11]. However, DWT has floating coefficient value that causing reversibility 
problem. Meanwhile, IWT is the next generation wavelet transformation which can solve the DWT problem 
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[12]. It is developed from the DWT using various lifting scheme to reserve the floating point within the range 
of integer value. 

This paper is proposed a new approach to embed the watermark on the integer wavelet domain using 
Chinese Remainder Theorem (CRT) algorithm. The low and high sub bands of IWT are used in this scheme. 
The high sub band is used to determine the embedding location on the low sub band and then the watermark 
bits are embedded using CRT algorithm based on remainders of mutual coprime integer c x and c 2 . The CRT 
is work on integer range that corresponds with the IWT and has lower requirement in space and complexity 
which can streamline the process. 


2. RESEARCH METHOD 

The proposed method is aimed to get better robustness within the same imperceptibility compared to 
the previous CRT method. The Haar IWT is used to decompose host image into four wavelet sub band of LL, 
HL, LH, and HH. The HH used to maps the embedding location on LL. The watermark bits are then 
embedded using CRT algorithm with c\ and c 2 value of 6 and 11 respectively. The robustness of the CRT and 
the proposed methods are measured under Gaussian filter, JPEG, and JPEG2000 compression using standard 
Normalized Correlation (NC). The next experiment is to measure imperceptibility the both algorithms using 
Structural Similarity (SSIM) by Wang et al. [13]. 

2.1. Haar IWT 

Haar IWT is an efficient integer wavelet that developed from Haar DWT using lifting scheme by Xu 
et al. [14 ]and many used in recent years by [15]-[17]. It is used simple round down operation to reserve the 
floating value in the decomposition process as shown in (1) until (4). Haar IWT is decomposed 2x2 
independence blocks of host image below 
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where I xy is an image pixel at row v and column y on each block. This process will convert the host image 
into four sub band with size of one quarter each. LL is a low frequency sub band that represents average of 
host image. The next sub bands HL and LH are the middle frequency that contains the horizontal and vertical 
feature. The last sub band HH has high frequency coefficient with the diagonal feature in it. 

The HH is used to determine the embedding location based on the absolute value of the coefficients. 
The highest coefficients value of the watermark bits numbers are then use to embed the watermark on the 
related position of LL. The reconstruction process is perform to build the watermarked image (/’) from the 
watermarked block. 
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In the reconstruction process the Haar IWT is added the integer value of 1 as a lifting scheme as 
shown in (5) until (8) to restore the rounded float value in the decomposition process. This is ensured perfect 
reversibility on the reconstruction process of the watermarked image. 


2.2. CRT Algorithm 

CRT is an algorithm that able to reduce the large size calculation within modulo operator to the each 
of modulus coprime factors [18]. It is able to simplify the determination of integer number from the 
remainder series using the coprime moduli [19]. Moreover, the CRT can be used as an embedding scheme 
through the modification of the coprime integers [20]. The 6-bit value of pixel p is defined as follows: 

pmodCj = r (9) 


where r, is the i- th remainder of p 
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The watermark bit b is embedded on the p through a modification process of 6-bit pixel value p 
based on the remainders r x and r 2 from the coprime integers c x and c 2 as shown in (9). The coprime integers 
are the mutual coprime values that use to generate the modulus m as shown in (10) until (12). The 
corresponding value a is determined such that the value of p is fulfilled. The modulus should be larger and 
closest to the maximum 6-bit value of 63 from the possible coprime integers. 
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Ho st Image Watermarked Image 



8 

Watermark Image 
Figure 1. Proposed method 


2.3. Embedding 

The p value is modified by the addition and subtraction of an incremental value and then the 
remainders are compared. It is done iteratively until they are fulfilled the condition of watermark bits. The 
detail embedding process is showed in Figure 1 as follows: 

a. Decompose the host image using (1) until (4) into wavelet sub band LL, HL, LH , and HH 

b. Read the number of watermark bits w 

c. Get the absolute coefficient value of HH and use the w-most significant value as pixel index k 

d. Store the 6-bit integer value as p and the 2 most significant bits value as q from the indexed pixels 

e. Calculate the remainders rl and r2 of indexed pixel using CRT: 


= p mod c t 

(13) 

r 2 = p mod c 2 

(14) 


f. Embed the watermark bit b according into p k based on this condition: 


f>r 2 ,b = 
\< r 2 , b = 0 
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g. Add the q into p to get the embedded sub band LL' 

LL — p® q (16) 

h. Reconstruct the LL', HL, LH, and HH to get the watermarked image. 


2.4. Extraction 

The k value is used for indexing the embedding location on the embedded sub band LL' in detail as 

follows: 

a. Decompose the watermarked image using (5) until (8) into wavelet sub band LL', HL, LH, and HH 

b. Get the embedded value p' using k from the indexed pixel of LL' 

c. Calculate the remainders rl and r2 of embedded value using (13) and (14) 

d. Determine the bit of watermark using the following condition: 


b = 


> >1 > r 2 
, r, < r 2 


(17) 


3. RESULT AND ANALYSIS 

In this paper, six standard images within size of 512x512 pixel in grayscale format and a binary 
image within size of 128x128 pixel are used as host images and watermarked image respectively as shows in 
Figure 2. 



(d) 


(e) 


(f) 



(g) 

Figure 2. The Host Images of (a)Baboon, (b)Boat, (c)Elaine, (d)House, (e)Jet, (f) Peppers, and 

Watermarked Image (g) 
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3.1. Robustness 

In this section, the robustness of the algorithms are measured under Gaussian filtering, JPEG, and 
JPEG2000 compressions with the standard sigma, quality factor, and compression ratio of 0.4, 90, and 3 
respectively. The result shows that the proposed method has outperformed the CRT in filtering and 
compressions test as presented in Table 1, Figure 3, 4, and 5. The CRT method result has the average NC 
value of 0.7129, 0.4846, and 0.6768 while the proposed method has higher NC value of 0.7902, 0.7473, and 
0.8163 in corresponding Gaussian filter, JPEG, and JPEG2000 compression test. This is due to the use of LL 
that has higher robustness than the direct embedding of CRT. It able to distribute the changes of Gaussian 
filter and both compression attacks to the down sampled of adjacent pixels evenly. The load partition process 
has successfully improved the robustness. 
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Robustness against Gaussian Filter 
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Figure 3. Robustness against Gaussian filter 


Robustuess against JPEG Compression 
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Figure 4. Robustness against JPEG compression 


Robustness against JPEG2000 Compression 
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Figure 5. Robustness against JPEG2000 compression 


3.2. Image Quality 

The next experiment is conducted to measure the images quality using standard SSIM due to the 
accordance with the human visual system. The Table 2 and Figure 6 show that the both algorithms method 
has similar performance in image quality measurement. The CRT and proposed method have average SSIM 





Robust Watermarking through Dual Band IWT and Chinese Remainder Theorem(Prajanto Wahyu Adi) 
















































568 n 


ISSN: 2302-9285 


value of 0.9979 and 0.9960 respectively. The proposed method is used the HH that contains the most 
significant coefficients which is directly proportional to the pixel frequency. 


Table 2. Images Quality 


Image 

CRT 

Proposed 

Baboon 

0.9994 

0.9993 

Boat 

0.9973 

0.9956 

Elaine 

0.9984 

0.9954 

House 

0.9980 

0.9943 

Jet 

0.9965 

0.9964 

Peppers 

0.9980 

0.9947 

Average 

0.9979 

0.9960 


Images Quality 



Figure 6. Images quality 


The location of the most significant coefficients are the mapped into LL to locate the coefficients 
that will be embedded. The selected coefficients are the highest frequency of down sampled pixel. They are 
able to reduce the degradation of structural similarity which is in turn maintains the imperceptibility. 


4. CONCLUSION 

This paper is proposed a new watermarking scheme to embed the watermark through dual band 
IWT. The watermark is embedded on the low frequency sub band which is indexed using the most significant 
coefficients of high frequency sub band. In robustness test, the CRT method result has the average NC value 
of 0.7129, 0.4846, and 0.6768 while the proposed method has higher NC value of 0.7902, 0.7473, and 0.8163 
in corresponding Gaussian filter, JPEG, and JPEG2000 compression test. Meanwhile the both CRT and 
proposed method have similar average SSIM value of 0.9979 and 0.9960 respectively in term of image 
quality. The result shows that the proposed method has higher robustness and able to maintain the image 
quality. 
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